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ABSTRACT 

Theories  of  rough  surface  scatter  and  gas  bubble  behavior  are  used 
with  the  Pierson-Moskovitz  wind-v:ave  spectrum  and  an  empirically- guided 
formulation  of  bubble  concentrations  at  sea  to  calculate  the  true  velocity 
gradient  and  losses  "at  the  surface".  These  values  are  then  entered  into 
Bucker's  wave  theory  solution  for  sound  propagation  and  leakage  in  a 
surface  duct.  Curves  of  propagation  loss  are  calculated  for  comparison 
with  ocean  test  data  obtained  with  the  SQS-26  sonar.  The  predictions  are 
shown  to  be  significantly  better  than  those  based  on  the  ernpirical  equations 
of  project  AMDS. 


I.  INTRODUCTION 

The  purpose  of  this  paper  is  to  show  that  physical  descriptions  of 
bubble  behavior  arid  sound  scatter  at  the  sea  surface,  extended  by  laboratory 
experiments  and  in-situ  measurements  of  near-surface  bubbles  have,  at  last, 
made  it  possible  to  give  reasonably  accurate  predictions  of  surface  duct 
propagation  from  a  knowledge  of  wind  speed  over  the  surface  and  duct  velocity 
profile • 

We  will  start  with  a  review  of  the  current  oceanographic  description 
of  the  sea  (Section  II),  and  surface  scattering  theory  (Section  III),  supple- 
menting the  basic  formalism  by  the  results  of  recent  laboratory  studies.  The 
theory  of  sound  dispersion,  scatter  and  absorption  by  bubbles,  simple  bubble 
kinetics,  the  extrapolated  results  of  in-situ  near-surface  bubble  measurements 
are  developed  in  Section  IV. 

(U)     Finally  (Section  V),  and  with  the  very  generous  cooperation  of  Dr. 
Homer  R.  Bucker  of  the  Naval  Research  Laboratory,  who  used  his  computerized 
wave  guide  solution  as  a  test  for  our  surface  loss  and  refraction  predictions, 
we  look  again  at  the  data  of  two  SQS-26  propagation  experiments  performed  by 
the  Key  West  Test  and  Evaluation  Detachment.  Bucker 's  wave  theory  solution  • 
is  modified  by  using  the  observed  wind  speed  at  the  time  of  the  experiment 
to  correct  the  velocity  gradient  for  bubble  presence  and  to  determine  the 
bubble  and  surface  losses  which  now  replace  an  empirical  one  db  surface  loss 
factor  that  he  used  to  obtain  excellent  data  fits. 


II.  THE  SEA 

We  assume  that  the  sea  surface  from  which  the  sound  is  scattered  is 
homogeneous,  and  that  the  temporal  variation  of  sea  heights,  C  ,   at  a  point 
is  Gaussian  distributed  with  mean  value  zero, 

"•■;■"  <o  -  o 
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and  with  vanance ,   o     , 


<C2>    -    a2 


The  surface  height  correlation  function  is 


c     „    A.     a     <C(0.0,0)  C(u.v.r))  (1) 

a2  a2 


where  u,  v,  T  are  the  spatial  and  temporal  lags. 

The  Fourier  transform  of  2  yields  the  three -dimensional  spectrum  of  the  sea 


X  (xx,x    .  0)  =  (2tt)'3  r  f        Z(u,v,r)  exp  [i(xu  +  xv  -nr)]dudvdr  (2) 

*T     U,V  X  y 


where       xv'x„   =    components  of  surface  wave  propagation  constant 
x     y 

O  =    angular  frequency  of  surface  wave 

We  shall  assume  the  gravity  wave  relation  between  angular  frequency,  propagation 
constant  and  acceleration  of  gravity,   g: 

O        =     xg  , 


This  permits  the  directional  spectrum  of  the  sea  to  be  specified  either 
in  the  two-dimensional  rectangular  coordinate  form 

X  =  *  (x,a)  6  (Vgx  -H)  (3a) 

or  the  polar  form, 

x  j  .♦  <XX<V    6(Vgx-n)  <*>) 

where  theDirac  delta  function, 6  (Vg*  -  Q)  ,   simply  selects  the  wave  that 
fulfills  the  gravity  wave  relation  between  frequency  and  propagation  constant, 
The  more  commonly  measured  frequency  spectrum   4>(Q)    is  now  obtained  by 
using  the  energy  density  coordinate  transformation 

/    >^KdKda  =    *(n)  dn 
0 

This  permits  the  frequency  spectrum  to  be  expressed  as  the  integral  of  the 
polar  form  of  the    X   spectrum  over  all  azimuth  angles, 

,3     .2f 
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*(fl)    =,—     \       *(xt  a)  da    .  .(U) 
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The  semi-empirical  form  of  the  frequency  spectrum  for  the  fully- 

Z 

developed  sea  that  is  most  widely  accepted  today      (Pierson-Moskowitz) 

is 


[-<(#] 


*(0)    -  ^      exp|-/5{-^J      |  .    (?) 


where 


a  =   8.1  x  10"3 

p  =    0.74 

"o  =g/W19.5         . 

W  =   Wind  speed  at  19. 5M  above  sea  surface 

g  ■    acceleration  of  gravity  in  consistent  units 


The  three  statistical  parameters  of  the  sea  surface  that  are  most  significant 
in  specifying  the  sound  scatter  are  the  mean  square  suface  height,  o-cf   the  mean 
square  surface  slope,  Z    $  and  the  surface  height  correlation  function,  C. 

For  a  locally-generated  sea  we  obtain  the  mean  square  height  by 
4  „+„~—  ~±i  n~   4-v.q    fVanv.jar>.rv   »no«*r:ni    fr.-r    fh~    fwl 1  v— dpvs»lCT>j=.H    r.?a. 


ttW4 


cr   =  C  0(0)  df2     -     2 
«L  4  Ba 


(6) 


For  the  mean  square  slope,  it  is  convenient  to  use  the  relation 

3 
determined  from  optical  measurements  of  the  sea 


E2    -    5.12  x  10"5W      +  0.003  +0.004  (7) 


where  W.,is  measured  in  cm/sec  at  a  height  of  kl   feet  over  the  surface. 
41 

Unfortunately  the  measured  speed  is  sensitive  to  the  height  of  the  anemometer j 

a 

however  corrections  can  be  made. 


Finally,  for  simplicity  and  because  of  the  general  lack  of  two- 
dimensional  sea  correlation  data,  we  will  assume  that  the  spatial  correlation  function 
is  of  Gaussian  form  and  isotropic: 


CO)  =  e"(S/L)  (8) 


where 

fl     y[T~         2 
X  =  vu  +  v 

and  •  L   =   surface  correlation  length. 

The  assumption  of  an  isotropic  Gaussian  function  makes  it  possible  to  calculate 
the  correlation  length  from  a  knowledge  of  rms  slope,  £  ,  and  rms  height,  q    , 
by  the  relation  that  is  valid  for  seas  of  lesser  roughness     : 


E  =  V2  s.  (?) 


III.     SOUND  SCATTERING  THEORY  AND  LABORATORY  STUDIES 


The  theoretical  solution  to  the  sea  surface  sound  scattering  problem 


usually  starts  with  the  Helmholtz  Integral: 

ikr„    .  .       /  ikr„v-i 


PQ 


4ff    ^S     L  F2        ^         P  bn     \    F2    / 


dS  (10) 


where    p_     =     the  scattered  acoustic  pressure  at  an  interior  field  point,  Q, 

of  an  enclosed  volume. 
S        =     integrating  surface  selected  so  that  acoustic  field  is  zero 

everywhere  except  on  the  ensonified  water  surface. 

k     •  =     2tt/\  =  acoustic  propagation  constant 

X        =     acoustic  wave  length  in  water 

r 2      =     distance  from  surface  scattering  region  to  field  position. 

n        =•     normal  to  surface  scattering  point. 
See  Figure  1. 
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After  applying  the  boundary  conditions  for  acoustic  pressure  and 

sous  tic  particle  velocity  at  the  sea  surface,  and  assuming  that  there  is  no 

ladowing  or  secondary  scattering  even  for  near-grazing  incidence,  we 

)tain  the  integral  for  the  scattered  pressure      (now  called  Ap2  s  P0)     from  a  given 

irface  of  scattering  area     A  A       in  terms  of  the   geometry: 


plane  AA 


iere     Ap9  =       scattered  pressure  for  a  plans  surface    A  A     over  which 

plane  there  is  negligible  curvature  of  the  incident  wave  front, 


ii        a  it  a      ik(r0  -  w.t)  /_ 

ikp.   AA  cos  0    e       2         i    /Ivi- 


AA  *     ,j*nsonif  ied  surface  area  over  which  p^  is  constant  and  the 

sagitta  of  the  wave  front  is  less  than  A/8  • 

p  at  /incident  sound  pressure  at  sea  surface. 

to.  ■     incident  sound  angular  frequency 

1  +  cos  9     cos  9     -  sin  8     sin  9?  cos  9 
cos  9     (cos  9    +  cos  9    ) 

8  =     angle  of  incidence  (measured  with  normal) 

9  =     angle  of  scatter  (measured  with  normal) 

9  =     azimuthal  angle  of  scatter  plane  with  respect  to 

incident  plane. 

Jfej  -     k  (ij  sin  9:  -  J_3  cos   91) 

k2  =     k  (JL    sin  $2  cos  9    +  i^  sin  92  sin  9    +  j^  cos  9   ) 

i 
•  i 

i.,  ,i0,i,,      are  the  unit  normals  in  the  x,  y  and  z  (depth)  directions, 
respectively. 

It  should  be  noted  that  the   solution -be ins  presented  assumes  that  the  illuminated 
surface  is  in  a  small  section  (JZX  <<  r±)   of  the  Fraunhofer  region  of  the  source 
field  and  that  the  field  position  is  in  the  Fraunhofer  region  of  the  illuminated 
area  (r2»VaA,   r2  >>  A  A/X  )•     This  point  is  emphasized  by  our  A  notation  for  the 
illuminated  area. 


The  mean  square  scattered  pressure,  which  in  general  includes  both 
coherent  and  incoherent  sound  intensity,  is  formulated: 


fcl^P^-Irt-^SSSS     <eiKz[C(x.y)-C(x'.y)]>ei[Kxu  +  Kyv] 

x  yxy 

(dxdx'dydy*) 
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where  K       -     JS1  -  Js2 

u  ■     X1  -  X 

v       »   y*  -  y 


Jpl  =  Ap2   ,         AP2*  , 

"  plane  plane 

*  =  complex  conjugate 

<   )  =  ensemble  average 


The  ensemble  average,   <    >  ,  can  be  easily  evaluated  for  our  surface 
with  its  two-dimensional  Gaussian  distribution: 

,    iK^CC-Ov  -Cka(cos  6    +  cos  6  )]2  [1  -C] 

\eZ  /    =    e  l  ^ 

Define  the  ACOUSTICAL  ROUGHNESS  OF  THE  SEA  SURFACE   =  g 

o 

g=    [ka(cos  6    +  cos  d2)~j*  (12) 

so  that  the  mean  scattered  intensity  can  be  written 


«vv>  ■  w^r-^SSSS  e^v+v>dxd,dydy. . 


(AA)' 

\  /         x  y  x'  y 


If  our  interests  were  to  be  restricted  to  low  roughness  (g  <  1)  at  this 

-g(l-C) 
point,  it  would  be  useful  to  use  the  approximation  e       ~  (e"^)(l+gC). 

This  would  permit  immediate  separation  into  a  coherent  term,  e~S,  and  an 

incoherent  term  which  the  integration  would  transform  into  the  two  dimensional 

surface  spectral  description,  eqn.  (  *fc)  ,  (Weiner-Khinchine  transformation). 


The  rarity  of  experimental  measurements  of  two-dimensional  surface  spectral 

descriptions  turns  us  away  from  this  approach. 

<>  t  Instead,   for  scatter  in  the  specular  direction,   it  is  convenient  to 

try  a  series  approximation  :to  the  assumed  isotropic  Gaussian  spatial  correlation 

5,6 
function  (Eq.   8).     The  integration         then  produces  the  relative  scattered 

intensity,    (compared  with  the  plane  surface  reflected  intensity)   specialized 

to  specular  scatter,   in  which      e=e        9     sn,Fsl,g  =  kA2  sin2<t, 

12        3 
where  <$>  =  9O-0^   =    grazing  angle. 

<Ap  Ap  *)  2        s  m 

plane  \        m=l  / 

The  two  terms  of  Eq. (13)  represent  coherent  scattered  sound  (the  e~s 
term)  and  incoherent  scattering  (the  second   term).  A  study  of  the  term 
in  parentheses,  defined  as 


s(g)         ,     •-*     2      £-  M> 

m=l 


shows  that  S  (g)    ->  g    for  g    «     1 

and  S(g)   ->  ~   for    g  »     1      . 

2 

Since  it  can  be  assumed  that  L    <  AA,  we  can  conclude  that,  in  the  limits, 

the  specular  scatter  for  a  Gaussian  sea  will  be: 

For  Small  Acoustical  Roughness,  g   <<   1:    (pp*)  =  e  (13>) 

ttL2 
For  Large  Acoustical  Roughness,  g  »    1:    (pp*)  =  "7~r~  - — 7~^ (16) 

2k  Z  AA  sin^# 

In  practice  the  coherent  torn  is  dominant  approximately  up  to  g~l>  and 

sometimes  beyond  that  roughness,  depending  or  the  beam  widt';.,  .Fig.  2 
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shows  the  relative  behavior-  as  a  function  of  acoustical  roughness. 
It  is  seen  that  for  small  acoustical  roughness,  the  intensity  of  specular 
scatter  is  the  factor  e~£  times  the  mirror  reflection  ,  and  the  concept 
of  a  reflection  coefficient  is  useful. 

However,  it  is  important  to  obssrve  also  that  the  value  of  dB  loss  per 
bounce  is  in  general  an  inadequate,  and  in  fact  misleading,  concept  for  the 
rough  surface  case.  When  g»  1  the  specific  value  of  the  scattering  area 
enters  into  the  asymptotic  value  of  (pp*)    which  will  therefore  change  from 
one  bounce  to  the  next  as  the  sound  wave  continues  to  diverge  and  scatter. 

When  the  Fraunhofer  conditions  are  violated  (large  beam  width)  it  is 
necessary  either  a)  to  add  the  scatter  contributions  from  the  several  areas 

A  A       or  b)  to  re-solve  the  problem,  assuming  a  diverging  wave  at  the 

7  ft 
very  beginning  of  the  derivation.  '   The  latter  attack  and  its  verification 

ft 

by  laboratory  model  show  that  for  specular  scatter  from  an 

exponentially  correlated  surface  there  is  virtually  no  difference  between  the 
Fraunhofer  and  Fresnel  solutions  if  g<l.  However,  for  g>l  the  scattered 
intensity  may  be  significantly  greater  than  predicted  by  the  single  Fraunhofer 
solution.  T.iis  is  because,  commonly,  significant  signal  will  exist  over  an  area 
containing  more  than  one  A  A  as  defined  in  our  derivation. 

Laboratory  models  have  been  used  in  order  to  supplement  the  capabilities 

9  10 

of  current  theory.  '    We  have  recently  studied  the  behavior  of  the  coherent 

component  for  cross-wind,  up  and  down  wind  specular  scatter  under  widely-varying 
laboratory  seas.  Our  "seas"  are  produced  by  combinations  of  a  paddle  and  five 
fans  blowing  air  over  the  surface  of  our  25'  x  6*  x  81   deep  anechoic  tank* 
In  this  way  we  obtain  water  waves  that  give  us  scaled  replicas  of  compound  sea 
and  swell  surfaces  that  are  realistically  time-varying,  anisotropic,  and  near- 
Gaussian.  The  probability  density  function,  spatial  and  te?  .oral  correlations 

*The  Gaussian  correlated  surface  would  give  similar  results. 


and  frequency  spectrum  compare  closely  to  ocean  measurements  ,  except  for 
scale  factors  for  frequency,  time  and  space  lag,  of  course. 

Figure  3  shows  three  of  the  laboratory  sea  and  swell  spectra.  The 
probability  density  functions  of  the  wave  height  for  the  three  cases  are  shown 
in  Fig.  hi   the  near-Gaussian  probability  density  function  is  very  close  to 
that  observed  at  sea.    Equation  (13)  states  that  the  ratio  of  the  intensity 
of  the  coherent  component  to  the  intensity  reflected  from  a  mirror-like 
surface  should  be  e"^.  Figure  5  shows  this  ratio  divided  by  e~&  over  four 
decades  of  acoustical  roughness j  for  g-  1  the  ordinate  is  essentially  unity. 
This  excellent  verification  of  theory,  when  tested  under  widely  varying 

laboratory  seas,  has  encouraged  us  to  apply  it  to  the  real  ocean 0  Other 

11  12 
laboratory  model  results  '  '  suggest  that  for  g <  1,  at  near  grazing  incidence, 

a  shadowing  correction  is  not  necessaryo 
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IV.   EFFECTS  OF  NEAR-SURFACE  BUBBLES 

Bubbles  near  the  sea  surface  can  modify  the  forward  scatter  of  sound 
in  two  ways:   (a)  By  near-surface  refraction;  the  presence  of  all  bubbles 
(resonant  or  not)  will  create  a  below-surface  medium  through  which  the  speed 
of  sound  differs  from  that  in  the  non-bubbly  water;  this  change  of  speed  will 
refract  incoming  rays,  modify  the  angle  of  incidence  at  the  rough-scattering 
surface  and  change  the  distance  betweeen  bounces  in  channel  propagation. 
(b)  By  resonant  absorption  and  scatter;  resonant  near-surface  bubbles 
can  cause  attenuation  that  is  generally  unresolvable  from  rough  surface 
scatter  loss. 

The  objectives  of  the  present  work  are  arbitrarily  restricted  to 
"surface"  losses  at   frequencies  below  16  kHz.   Nevertheless,  because  bubble 
refraction  depends  on  the  total  volume  of  air  in  bubble  form,  it  is  necessary 
to  consider  the  numbers  of  near-surface  bubbles  of  ail  radii. 

In  SDite  of  the  great  interest  in  the  subject  for  almost  30  years, 
the  literature  of  ocean  bubble  data  is  very  thin:  In-situ  measurements  of 
bubbles  in  off-shore  breaking  waves  have  been  reported1-*;  bubble  generation 
and  persistance  by  an  artifical  wind  over  a  laboratory  tank  filled  with  sea 
water  has  been  described  ;  inferred  bubble  densities  have  been  calculated 
from  certain  AMOS  data  ;  in-situ  acoustical  measurements  in  60  ft  of  iso- 
thermal water  in  Mission  Bay,  San  Diego  have  been  summarized  ;   a  laboratory 

study  of  surface  bubbles  has  yielded  data  produced  by  colliding  waves  of 

17 

sea  water e    It  will  be  necessary  for  us  to  extrapolate  beyond  these  bare 

measurements  available. 

A.  Dependence  on  Radius 

Because  of  experimental  difficulties,  in  no  case  ar2  there  published 
in-situ  measurements  of  ocean  bubbles  of  radius  greater  than  200 

13 


microns  (resonant  frequency  less  than  approximately  15  kHz)  or  less  than 
roughly  20  microns  (resonant  frequency  greater  than  150  kHz) .   Most  of  the 
available  data  fit  approximately  the  solid  line  in  Fig.   6  >  plotted  on 
the  basis  of  fractional  volume  of  air  in  bubble  form  to  water  u(R)dR.   From 
the  fractional  volume  it  is  easy  to  calculate  the  number  of  bubbles  per  unit 
volume,  in  radius  increment  dR,  by  using 


>(D  dR    -   a«Jtt  (1?) 

4/3nR 


There  is  some  evidence  that  the  bubble  populations  break  up  roughly  into  two 

groups:   those  smaller  than  approximately  60  micron  radius  postulated  as 

attached  to  motes  and  due  principally  to  biological  activity,  cosmic  rays, 

captured  aerosols  and  photosynthesis  cf  phytoplanktoiij  and  the  larger  bubbles 

which  we  assume  come  from  breaking  waves  (and  decaying  matter  on  the  sea 

floor  in  shallow  water) . 

,v-    To  calculate  the  near-surface  refraction  due  to  bubbles  we  have  had  to 
extrapolate  (dashed  lines)  the  data-based  curves  of  Fig.  6,  which  were 

measured  with  incipient  breaking  waves  at  wind  speeds  U-6  knots  at  3.3  meter 

depth.  For  radii  smaller  than  30y  we  extrapolate  to  one  micron  radius  at  the 

-1/2 

h^B.     '        dependence  that  has  been  found  for  the  smallest  measured  ocean 

bubbles   and  in  tap  water  j  below  1  micron,  at  slope  K^R  following  behavior 

19  20 
loosely  suggested  by  laboratory  cavitation  studies.   *    For  bubbles  larger  than 

60n  we  have  used  the  KiR   dependence  observed   to  100  microns  and  have 


1U 


extrapolated  to  largest  radii  at  Blanchard  and  Woodcock's  slope  of  K  R 
The  radius  at  which  one  shifts  to  the  suggested  slope  for  largest  bubbles 
is,  admittedly,  somewhat  arbitrary.   That  slope  starts  at  about  50  microns 

in  the  breaking  wave  data;  however,  a  re-analysis  of  Glotov's  data  shows  that 

-3 

the  K  R  "  behavior  is  confirmed  for  R  >  120p  .   The  constants  for  the  slopes 

of  Fig.  6  are  shown  in  Table  1. 


TABLE  ONE 

Estimated  Fractional  Bubble  Air  to  Water  at  3.3M  depth  for  Winds  Approximately 
I  6  Knots 


Radius  Range        Dependence        Constant  [dimensions]    /  u(R  )  dR 
Microns  u(Rn,3.3M) 


RQ  <_   1.0  K^0  X  -  1.5  x  10"9  [y-1]  1.5  ::  10"9 

1*°  iR0£  30-  K2R0~1/2  K2  =  1*5  X  10"9  ^~1/21  13'8  x  10~9 

30.   <R  <  60.  K3R0_1  K3  =  8*5  X  10"9  ^°3  5*8  X  10~9 

60.   :  RA  :  100.  K.R."*"1  K,  =  2.5  x  10~12  [y"2]  7.7  x  10"9 


-  0  -    *         4  0        ^4 
i  R0  K5Ro'3       K5 


100.  I  Rn  K_R  ~3       K_  -  2.5  x  10~4    [u+2]    12.4  x  10  9 


—8 
Total  Effective  U  =  4.1  x  10 
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B.   Dependence  on  Depth 

For  our  purposes,  interest  lies  in  the  bubble  densities  (and  fractions) 
as  a  function  of  depth  at  a  fixed  frequency,  rather  than  for  fixed  radius. 
Studies  in  well-mixed  ocean  water  suggest  that  for  bubbles  larger 

than  approximately  30  microns  (sea  level  resonant  frequency  less  than  100  kHz) 
buoyant  forces  and  turbulent  entrainment  determine  the  bubble  life  time 
rather  than  gas  diffusion  through  the  bubble  walls.   If  we  formally  make  the 
simplifying  assumption  that  these  bubbles  can  be  lost  only  at  the  ocean  surface, 
then  bubbles  of  surface  (3.3  Meter  depth)  resonant  frequency  f-,  radius  R^, 

will  be  found  at  depth  Z  in  the  same  numbers,  but  at  the  isothermally* 

1/3 
compressed  radius,   R.   =  R.3     and  resonant  frequency  f    where 

p    -    1.33/(1  ♦  0.1Z). 

For  resonant  air  bubbles   kRn  =  ^   f&?    where  c  «  LSxlQ5  cm/see 

y-  1.4,  P  =  106  ^WlOdyne/cm2^  m  1.0  gm/cm3,  Z  =  depth  in 

meters .  


The  bubble  fraction  at  depth  Z  will  be 

u(RQz,Z)   =  n(R02,Z)(4/3  .R^3) 

Analytically,  our  assumption  for  large  bubbles  is:  from  the  surface 
to  depths  where  the  entrained  bubble  drag  force  is  equal  to  or  greater  than 
the  buoyant  force 

n(RQ,3.3M)  dR  =   n(RQZ,Z)  dR  , 
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80  that  following  isothermal  compression, 

9  --r  — 

u(RQz,Z)   =   [u(R0,3.3M)][B].  (18) 

These  bubbles,  however,  have  a  resonant  frequency  which  is  not  f  ,  but, 
because  of  changed  radius  and  pressure,  is 

f    -   f   B~5/6 

foz        fo  e 


We  can  now  calculate  the  gas  fraction  of  bubbles  at  the  desired 
resonant  frequency  f_  at  depth  Z,   u(f  ,Z)  ,  by  assuming  that  the  frequency 
dependence  within  a  given  band  of  bubble  sizes  is  the  same  as  at  the  surface, 
that  is, 


u(f  ,3.3M)   =  k  .  f_      and     u(f„_.Z)   -   k  „  f  „ 
0  n'-J     u  W        nZ  Ct. 


where  n  is  the  frequency  dependence  within  a  particular  frequency  band. 
Therefore 


u(f0'Z)  = 

n        -,5/6," 


u(foz,Z)       0  01 


(fn/f   )        -     (r1  )       •  (19) 


Finally,  from  (  18)  and  (19) 


u(f0,Z)   =   (B5/6  n+1)  u(fQ,3.3M) 


and,  taking  account  of  the  change  in  radius, 


(20) 


n(fo,z)   .   (B5/6n+5/2)  n(f0,3.3M)  (21) 
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The  predicted  depth  dependences  for  the  three  most  important  bands 
of  surface  radii  are  shown  in  Table  2# 

TABLE  2 

Fractional  Air  in  Bubble  Form  and  Number  of  Bubbles  per  Unit  Volume 
as  a  Function  of  Depth  for  Large  Bubbles 


u(f0,Z)/u(f0,3.3M)   N(fQ,Z)/N(f0,3.3M) 

fi5/6  n+1  .5/6  n  +  5/2 

■  p  ■  p 


+1  6 1V6  B10/3 

-1  81/6  65/3 

>  100  less  than  32   +3  B  B5 


Ro 

fo 

microns 

kHz 

30  to  60 

60  to  100 

60  to  100 

32  to  60 

It  is  clear  that  the  assumptions  that  lead  to  the  values  in  Table  2 
can  only  be  considered  as  a  perturbation  analysis  of  the  change  of  bubble 
populations  with  depth  since  they  imply  inconsistencies,  particularly  at  the 
edges  of  the  bands,  as  depth  is  increased.   Nevertheless,  the  only  bubble 

density  data  available         plotted  against  the  absolute  pressure,  1.333   , 

+2 
for  four  frequencies  between  40  and  96  kHz  show  slopes  of  approximately  3   , 

1.7       3.3 

in  reasonable  agreement  with  the  predictions  of   8     to   8 

When  we  plot  the  frequency  band  contributions  to  the  total  bubble  fraction, 
(  /  u  R  dr) ,  as  a  function  of  total  pressure  parameter   8   ,  the  total  bubble 

fraction  U(z)  at  wind  speed  6  knots  is  found  to  vary  (with  less  than  5%   deviation) 

+2  +1 

as   8    from  the  surface  to  approximately  10M  depth,  blending  into   6 

at  greater  depths  (near  20  Meters).   Therefore  the  dpeth  dependence  that  we 
will  use  for  the  near-surface  refraction  calculation  is: 

U(Z)  =  4.1  x  10'8  8+2  (22) 
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C.   Dependence  on  Wind  Speed 

It  has  been  found    that  the  number  of  bubbles  shows  "an  abrupt 

increase  ...  approximately  ...  exponential"  for  winds  above  about  6  knots. 

21 

The  comment  agrees  roughly  with  the  curve  of  foam  coverage  vs.  wind  speed 

which  can  be  written  approximately  as: 

U(Z=0,W)   =  U(Z=0,6  knots)  [1  +  10~2(W/10)2]  (23) 

where  W  =  wind  speed  in  knots. 

There  is  a  second  wind  effect:   the  larger,  wind-generated  bubbles  will  be 
brought  to  depths  by  increased  turbulent  entrainment*  .  A  complete  knowledge  of 

this  behavior  may  be  possible  when  there  is  a  better  understanding  of  under- 

22 
water  motion  as  a  function  of  wind  speed,  insolation  and  depth. 

However,  the  factor  of  bubble  entrainment  by  water  motion  can  be  crudely 

evaluated  by  the  following  simple  argument: 

A  bubble  will  be  drawn  down  from  the  surface  when  the  drag  force  is 

greater  than  the  buoyant  force. 

The  coxivective  velocity  in  the  medium  will  probably  be  Gaussian  distributed. 

Nevertheless,  the  mean  downward  convective  velocity  will  increase  with  wind 

speed  as  does  the  velocity  spectrum   and  the  depth  of  the  well-mixed  layer. 

We  assume  that  the  mean  convective  velocity  is  proportional  to  the  wind 

speed  and  that  the  mean  number  of  bubbles  drawn  downward  under  a  given 

condition  of  underwater  turbulence  is  proportional  to  the  mean  velocity;  the 

dependence  on  radius  will  be  ignored  at  this  time.  This  "zero  order" 

estimation  of  the  effect  of  turbulent  entrainment,  combined  with  eq.(23) 

for  the  increased  foam  available  leads  to  the  wind  dependence: 


U(Z,R  ,W) 

=(w/6)  [l  +  io  Wio)2] 

U(Z,R  ,6  knots) 


1$ 


(2)4) 


D.  Summary  of  Bubble  Dependences 

Combining  the  dependences  on  radius  of  bubble,  bubble  depth  and  wind 
speed,  we  get  the  total  bubble  fraction: 

U(Z,W)  »  [U.l  x  10~8£2]  [W/6]   [1  ♦  10"2(WA0)2]  (25) 

It  is  appropriate  to  point  out  that  other  bubble  factors  exist  in 
addition  to  the  effects  that  we  have  considered.  Two  significant  elements 
are  the  dependence  on  daylight  and  precipitation;  others  are  dependence  on 
water  temperature  and  gradient,  percentage  ga3  saturation,  marine  biological 
activity.  As  examples:  Blanchard  and  Woodcock  have  shown  that  snow  and 
rain  produce  very  large  numbers  of  bubbles.  Their  penetration  in  laboratory 
experiments  was  no  more  than  2  cm  below  the  water  surface  but  it  would  be 
expected  that  these  bubbles  would  penetrate  much  farther  in  the  presence  of 
convection  cells  accompanying  storms  at  sea.  Small  bubbles  (less  than  about 
60  v)   are  about  twice  as  populous  in  the  daytin.e  as  at  night.  It  has  been 
postulated  that  the  increased  numbers  may  be  due  to  photosynthetic  action. 
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E.      Theory  of  Dispersion  and  Absorption  by  Bubbles 

Following  classical  developments  ^5      an  incident  acoustic  pressure , 
of  wave  length  very  large  compared  to  the  bubble  radius  R,  having  a  velocity 
potential  $  e        ,  will  produce  a  total  field  , 


N 


0(ir)    =    0.  (ir)    +   1       (-f-) 
1  n  =  1    x  r  ' 


ik  I  ir  -  ir   | 
e  n 


(26) 


where    ir 


n 


n 


R 


0) 


o>. 


coordinate  position  of  scattering  region 

th" 


distance  between  n      bubble  and  center  of  scattering  region 
2tt/X  =  propagation  constant  in  bubble-free  water 
Rtf) 


n 


(w0/w)     -  1  -  i6 

—  bubble  radius 

-  incident  sound  frequency 

=     resonant  frequency  of  bubble 
=     damping  constant  of  bubble 


The  average  over  all  possible  sizes  and  configurations  of  bubbles  is 
<*<«»  -  #1  OD    +/J/G  Orj    <0(«rn)  >     ^JH^    dVn 


(27) 


where 


G  00    -    f    In  (ir.R)3 
n  R 


R  d  R 


Suo  (R)/«3     -  1  -i6M 
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(u) 


Equation  (27)is  the  integral  equation  formulation  of  the  differential 


equation 


V2    <0>    +    k2   <</)>    =     -4ttG<0>         or 


V2    <tf>>    +   kb2  (r)    <0>    =     0 


(28) 


where 


kb(r)    =    k 


b 


AttGI 

k2J 


«  '  k  + 


277       r 


n  (r,R)    R  d  R 


R       K/w)"  -  1  -i6(R) 


kb(r)  47]>G 

The  quantity,  — : =  1  +  — ~  ,  is  the  index  of  refraction  in  the  bubbly 

region.    Because  of  the  presence  of  the  damping  term,  the  index  will  have 

an  imaginary  part  which  comprises  the  attenuation  due  to  absorption  and 

scatter  in  the  bubbly  water  in  addition  to  the  real  part,  defining  the  local 

speed  of  sound.    These  terms  are: 


Re:     {k    (r)}     =     k  + 


2£       f 
k       J ' 


■  »4fe- 
[(^)2  -  >] 


- 1 


d  R 


+    6' 


(29) 


K 


6  R   n  (R)  d  R 


[(» 2  - .] 


+     6' 


(30) 


22 


F#  Predictions  of  Refraction  by  Bubbles 

The  dispersion  relation  Eq.   (29) »  since  it  predicts  a  decrease 
of  sound  speed  caused  by  bubbles  of  resonant  frequency  greater  than  the 
incident  frequency  and  an  increase  of  speed  for  bubbles  with  lower  resonant 
frequencies,  would  require  precise  information  of  bubble  populations  for 
an  evaluation  for  all  sound  frequencies.   However,  if  our  interest  is  restricted 
to  sound  frequencies  less  than  16  kHz,  integration  of  Fig.   6  shows  that 
bubbles  under  that  resonant  radius  (R~  <  200.  microns)  represent 
approximately  93%  of  the  total  fractional  air  volume.   For  the  calculation 
of  refraction  at  frequencies  less  than  16  kHz  we  consider  only  the  effect  of 


bubbles  of  resonant  frequency  greater  than  the  driving  frequency. 

2 
When  w0  »  u  ,  we  drop  6  in  the  denominator  of  Eq.  (29l,i 

w0/w  =  R0/R  and  the  definition  of  u  (R)  to  obtain  the  simplified  form: 


Re:   {kb}   =  k  (l.+ 


3  U   \  ,  . 

w  «  «n       (3D 


2k2R02; 


Since  the  expression  in  parentheses  is  nearly  unity,  we  find 

c,   =  c   1  -  8.0  x  103  U   (1  +  Z/10)"1  (32) 


1 


where    U(Z,W)  is   given  by  equation  (2£). 
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The  significance  of  this  departure  of  the  speed  from  the  bubble-free 
value  is  evident  when  we  insert   Equation  (25)       and  look  at  the  gradient 

of  the  near-surface  propagation  speed.   For  the  reference  wind  speed  of 

-4 
6  knots  the  gradient  of  the  velocity,  proportioanl  to  (1  +  0.1Z)    f  has  a 

surface  value  of  0.26  sec   ,  decreasing  to  0.016  sec  '"  by  10  Meter  depth. 

Since  the  velocity  gradient  due  to  pressure  in  isothermal  water  is  0.018  sec"" 

our  conclusion  must  be  that  ray  curvature  in  approximately  the  top  10  meters 

is  very  strongly  influenced  by  bubbles  even  for  winds  as  low  as  6  knots. 

Because  of  bubbles,  rays  approaching  the  surface  will  refract  upward, 

meeting  the  surface  in  a  cusplike  ray  path,  and  making  more  bounces  for  a 

given  horizontal  range  than  for  a  bubble-free  region. 


2k 


\  -) 


G.     Predictions  of  Attenuation  due  to  Bubbles 

The  sound  pressure  attenuation  formula  for  the  bubbly  water  can  be 
written  as 


^-'-Ofx     i(wt-k,  x) 
p     =     p    *  e  e  b 


where 


a      =     Im:    {k :    } 


Using  Eq.(30)and  rewriting  in  decibel  notation,  we  get  the  attenuation 


a,      «     8.68  a 

D 


/*    (6/6.)  R      n(R,Z,W)   dR 


(8.68)  (2ir)J 


CO 


0 

(a) 


-1 


+  6* 


(33) 


where         6       =     radiation  damping  constant   =    k  R 


ftt 


1.36  x  10"2   ^1+Z/lO     for  air  bubbl 


es 


26 
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Integration         yields 


a       «=      [8.68ir2R  3n(Ro,Z,W)]/6r 


which  can  be  re-expressed  as 

,9  „L„,.,nx-l/2 


av  =  1.50  x  10'  (l+Z/10)'-"-^  u(fn,Z,W)   db/Meter       (3k) 
b  U 


The  conclusion  of  the  calculation  of  "surface"  loss  in  db/bounce 
due  to  bubbles  requirer  integration  of  Eq.  (3h  )  over  the  surface  duct  ray 
path  of  varying  curvature  and  varying  attenuation: 

0 

*b  "  2    I  %   dS 

where  L  is  the  depth  at  which  the  ray  direction  is  horizontal.   Inserting 
a,   and 


Sln9  2j2«.(L-2)     2  V  2(L-Z)G 

where   $.  =  local  radius  of  curvature,  C-(Z)  is  the  total  local  gradient,  and 
c  =  local  velocity,  M/sec,  we  obtain 

o   , 0    u(f   Z,W)         dZ 

A^  =  1.06  x  10y  FT   /    ^-ry^ U2      . /2 (35) 

*  4     L     (L-Z)1/2  (1  +  0.1Z)1/2  G1/2(Z)   ' 


For  frequencies  below  16  kHz,  the  expression  for  the  bubble  fraction, 
from  Table  2  and  Section  C,  is 

«(*0.Z*W)  *  0.8  x  10"2U  fQ3  ^7/2W[l  +  10-2(W/10)2J  (36) 


where  W  is  in  knots  and  fQ  is  in  Hz0 
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V.  SURFACE  DUCT  PROPAGATION  -  COMPARISON  WITH  EXPERIMENT 

1 
Bucker  has  recently  solved  the  surface  duct  propagation  problem  by 

a  wave  theory  approach.  His  nodal  components,  which  correspond  to  up  and 

down  propagating  waves,  with  appropriate  ray  directions,  lose  energy  when 

they  interact  with  the  surface  and  by  leakage  through  the  thermocline.  Our 

proposal  here  is  to  introduce  the  results  of  our  development  in  Sections  III 

and  IV  into  Bucker' s  solution.  Two  changes  are  necessary:  the  velocity 

gradient  in  the  duct  will  now  be  due  to  bubbles  as  well  as  temperature  and 

pressure  effects ;   the  loss  at  the  surface  will  he   made  up  of  the  losses  per 

bounce  due  to  below-surface  extinction  of  sound  by  bubbles  and  simple,  at-surface, 

coherent  specular  scatter  (because,  in  general,  g<l),  calculated  for  each 

contributing  mode.  The  surface  "reflection  coefficient"  is  thereby  directly 

attributable  to  the  wind  velocity  which  is  the  sole  source  of  the  rms  wave 

height  in  the  examples  chosen  and  a  principal  cause  of  the  bubbles. 

First,  let  us  look  at  how  the  losses  vary  with  wind  speed  and  ray  angle 

of  incidence.  Our  loss  calculations  are  based  on  eqs.  l£,  3$ >   36.  Our  ray 

paths  are  calculated  from  the  bubble-caused  gradient  using  eqs.  2£,  32  together 

with  a  temperature  plus  pressure-caused  gradient  of  0„0U8  sec  .  Figure  7 

illustrates  the  relative  contributions  of  surface  loss  and  integrated  bubble 

loss  for  a  range  of  winds  with  a  ray  that  is  horizontal  at  UO  M,  ^  =  U.38°. 

below 
The  frequency  isA  16  kHz.  It  is  observed  that  the  bubble  contribution 

dominates  only  for  low  wind  speeds. 

*In  the  data  report   the  "seaman's  eye"  estimates  of  the  peak-to-peak  wave 
height,  H,  are  in  good  (8C>)  arrcement  with  the  Calculation  from  equation  (6) 
when  it  is  assumed  that  H  =  6o~  ,  for  Gaussian  wave  systems.  However,  our 
laboratory  experiments  suggest  that,  in  general,  the  total  height  due  to 
swell  and  "sea"  should  be  used  to  determine  cr. 
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Figure  8  shows  that  the  bubble  loss  per  bounce  i3  not  very  sensitive 
bo  the  depth  of  the  horizontal  ray  (or  angle  of  incidence  at  the  surface); 
the  loss  per  bounce  at  the  surface,  however,  increases  as  the  square  of  the 
sine  of  the  grazing  angle. 

For  higher  frequencies  it  would  be  expected  that  the  neglect  of  the 
Incoherent  surface  scattering  term  would  no  longer  be  possible;  the  forward 
scatter  loss  would  reach  an  asymptotic  value  which  would  be  smaller  than  the  increased 
nubble  loss.  For  frequencies  lower  than  our  example  the  coherent  scatter 
approximation,  e"^,  becomes  even  better,  and  would  be  applicable  to  higher 
sea  states;  the  bubble  loss  would  be  expected  to  continue  to  be  smaller  than 
the  surface  loss. 

The  experimental  propagation  data  for  our  study  are  again  from  Runs  73 
ind  99   of  the  operational  tests  of  the  SQS-26  sonar  '  in  which  both  source 
and  receiver  were  in  the  duct.  The  wave  heights  are  such  that,  g<l  for  ail 
affective  modes  of  the  duct  solution. 

Since  Bucker's  duct  solution  assumes  constant  velocity  gradients,  we 
lave  had  to  approximate  the  calculated  total  gradient,  which  goes  as  p    , 
sy  splitting  it  into  constant  sections  G-^  and  Gp  (.Fig.  9).  The  near-surface 
section,  G]_  ,  has  been  given  a  value  twice  the  uncorrected  gradient  due  to 
temperature  and  pressure  and  the  remainder,  Gp  ,  has  been  assumed  to  have  no 
bubble  influence.  The  values  are  given  in  Table  3.  In  this  way,  the  duct 
propagation  loss  curves  have  been  recalculated  by  Bucker.  They  compare  well 
dth  the  data  (Figs.  10,  ll).   The  dashGd  lines  in  Figures  10  and  11  are  the 
much  poorer  predictions  based  on  the  crruations  of  rroject  N-TDZ     / 
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TABLE  3 
Data  for  Two  Surface  Duct  Experiments 


Run  No. 

vs 

Zt 

Gl 

G2 

W 

M/sec 

M 

sec"^- 

sec"l 

knots 

73 

1513.1 

36.6 

0.096 

0.0^8 

17 

99 

1510.9 

51.8 

0.102 

0.031; 

15 

where   V  =  velocity  at  surface 

Z*  ■  depth  of  thermocline 

G-^  *  velocity  gradient,  surface  to  8  yards 

Gg  ■  velocity  gradient,  8  yards  to  thermocline 


Other,  higher  values  of  the  assumed  bilinear  gradient  and  a  closer- 
fitting  tri-linear  gradient  have  also  been  considered.  However,  the  resulting 
propagation  prediction  curves  show  excessive  losses  compared  to  experiment. 
Since,  at  these  wind  speeds,  the  bubble  loss  is  less  than  the  surface  loss, 
the  velocity  gradient  can  be  a  major  influence.  This  suggests  that  our  total 
bubble  populations,  U(fQ,Z,W),  may  be  too  high,  perhaps  by  a  factor  of  two, 
possibly  due  to  the  wrong  wind  speed  dependence  having  been  assumed. 
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VI.  CONCLUSION 

The  purpose  of  this  paper  has  been  to  show  that  for  a  large  range  of 
ocean  surface  roughness,  g < 1,  it  is  possible  to  predict  surface  duct 
propagation  by  using  the  simplest  surface  scattering  theory  (coherent  term 
only)  and  estimates  of  bubble  presence  as  inputs  to  the  Bucker  modal  solution. 

Nevertheless,  it  is  the  necessity  of  knowing  more  about  the  ocean 
parameters  of  sound  propagation  that  has  been  demonstrated  by  this  paper. 
It  is  hoped  that  the  present  work  will  serve  as  a  spur  for  more  fundamental 
studies  to  replace  the  speculative  assumptions  that  we  have  been  forced  to 
make  about  bubble  dependence  on  wind  speed. 
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